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Abstract 

We consider the lattice dynamics in the half-space. The initial data are random 
according to a probability measure which enforces slow spatial variation on the linear 
scale . We establish two time regimes. For times of order e~'^ , < 7 < 1 , locally 
the measure converges to a Gaussian measure which is time stationary with a covariance 
inherited from the initial measure (non-Gaussian, in general). For times of order , 
this covariance changes in time and is governed by a semiclassical transport equation. 
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1 Introduction 



The paper concerns a mathematical problem of foundations of statistical physics and contin- 
ues the work [9] devoted to the derivation of a limiting "hydrodynamic" (Euler type) equation 
from the Hamilton dynamics. We refer the reader to [2, 4, 14, 15] for a detailed discussion 
of the results and methods on this problem. 

As the model we consider the harmonic crystals in the half-space Z'l = {z e : zi > 0} . 
In the harmonic approximation, the crystal is characterized by the displacements u{z,t) e 
R" , z G Z'l , of the crystal atoms from their equilibrium positions. The field u{z, t) is 
governed by a discrete wave equation. 

The derivation of hydrodynamic equations is connected with the problem of convergence 
to an equilibrium measure. Hence, the first step in our inverstigation is the proof of such 
convergence. This step was done in [10]. We assume that a probabihty measure giving 
the distribution of initial data has some mixing properties. If denotes the time-evolved 
measure at time t , then the limit 

lim lit = /Xoo (1-1) 



t— »oo 



is established, where yUoo is an equilibrium Gaussian measure. (The precise formulation of 
this assertion is given by Theorem 2.5). In [5, 6], we have analyzed the long-time convergence 
to an equilibrium distribution for systems described by partial differential equations in R*^ . 
In [7]-[9], we extended the results to harmonic crystals. In the above-mentioned papers the 
systems were considered in the entire space. In [10], the dynamics of the harmonic crystals 
is studied first in the half-space . 

To derive the hydrodynamic equation we apply the special so-called hydrodynamic limit 
procedure. Given a matrix function {i?(r, •),r e R'^} (so-called "spectral density matrix 
function" in the terms of R.L. Dobrushin and others, [3]) we consider a family of measures 
{//q,£ > 0} which satisfies the following conditions: (i) For any r e R^ , the covariance 
Qs{z,z') of the measure /j^ at points z,z' £ close to the [r/e] is approximately (as 
e — > ) described by i?(r, ■) ; (ii) the covariance (5e(-2, z') vanishes as \z — z'\^oo uniformly 
in s (see conditions VI and V2 in Section 2.2 below). Given nonzero r e R and r G H'^. , 
we study the distribution /^^/^^ of the random solution u{z,t) at time moments r/e and 
close to the spatial point [r/e]. We establish the limit 

}yr/e,r = f'?,r, (1-2) 

where /x^^ is a Gaussian measure (see Theorem 2.14). In particular, we derive the explicit 
formulas for covariance matrix q^,.i^ — z') of the limit measure /i^j. . These formulas allow 
us to conclude that the matrix function q^^y.{6) evolves according to the following equation: 

drU,M =iC{9)Vn{9) ■ VrU,M, reZX, r > 0, (1.3) 

where C{9) = ^ ^(^Q'j ^ ^ ' ^^oughly, ^l{6) is the dispersion relation of the 

harmonic crystal. The boundary and initial conditions for (1.3) are written in terms of the 
function R{r, ■) . The equation of type (1.3) is called a hydrodynamic (or Euler) equation 
(see [3]). The result (1.3) is a continuation of the works [3] and [9]. In [3], the problem has 
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been studied for the infinite chain of harmonic oscillators on one-dimensional lattice . In 
[9], the result has been extended to the many-dimensional case. 

In phonon physics it is standard practice to use the Wigner function W{t, r, 9) as density 
of phonons with wave number 9 at location r and at specified time t . W evolves according 
to the semiclassical energy transport equation 

dtW(t, r, 9) = -VVL{9) ■ VrW{t, r, 9) (1.4) 

(see Theorem 2.12 and Corollary 2.13). W{t,r,9) at fixed r,t are expressed by the co- 
variance q^r{9) which is invariant under the lattice dynamics. Thus (1.3) or (1.4) can be 
understood as the equations governing the motion of the parameters which characterize the 
locally stationary measures. 
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1.1 Model 

We study the dynamics of the harmonic crystals in , d > 1 , 

u(z, t) = - J2 (^(^ - - ^(^ - ^0) t), ^ e Z^, teR, (1.5) 

with zero boundary condition, 

u{z,t)U,=o = 0: (1.6) 

and with the initial data 

u{z, 0) = Uo{z), u{z, 0) =ui{z), z e (1.7) 

Here Z^ = {2; G Z'^ : 2:1 > 0} , z = (—2:1, 2:2, ... , Zii) , V(^z) is the interaction (or force) ma- 
trix, iyu{z)) , kj = l,...,n, u{z, t) = {ui{z, t), . . . , ujz, t)) , Uo{z) = iuoi{z), Uoniz)) e 



R" , and correspondingly for Ui{z) . To coordinate the boundary and initial conditions, we 
assume that uo{z) — Ui{z) — for 2:1 = . 

Write Y{t) = {Y'{t),Y\t)) = {u{- , t) , u{- , t)) and Yo = (Y.'X) = {M-):M-)) ■ Then 
(1.5)-(1.7) becomes the evolution equation 

Y{t) = A+Y{t), teR, zeZi, Y'{t)l,=o = 0, Y{0) = Yq. (1.8) 

Here A+ ^ I I] with V+uU) E {V(z - z') - Viz - ~z'))uiz') . 

Let us assume that 

V{z)^V{z), where z = (-^1,5), z ^ {z^, . . . , za) e Z''-^ (1.9) 

Then the solution to the problem (1.8) can be represented as the restriction of the solution 
to the Cauchy problem with odd initial data on the half-space. More precisely, consider the 
following Cauchy problem for the harmonic crystal in the entire space Z*^ : 

yi^z.t)^- E V{z- z')v{z',t), z e Z'^, ieR, 

2'ezrf (1.10) 
v{z, 0) = Vo{z), v{z, 0) =^1(2;), z e Z"^. 
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Write X{t) = {X\t),X\t)) = {v{;t),v{;t)) and Xo = {X^,X^) = {vo{-),v,{-)) . Then 
(1.10) becomes 

X(t)^AX(t), teR, X{0)^Xo. (1.11) 



Here A = 




, where V is a convolution operator with the matrix kernel V . 



Assume that the initial data Xo{z) form an odd function with respect to Zi & , i.e., 

let Xq(z) = ~Xq{z) . Then the solution v{z,t) of (1.10) is also an odd function with respect 
to ^1 e . Restrict the solution v{z,t) to the domain Z^|_ and set u{z,t) = v{z,t)\zi>o ■ 
Then u{z,t) is the solution to the problem (1.5) with the initial data Yq^z) = Xq{z)\z-i^>o . 

Assume that the initial data Yq for (1.8) belong to the phase space 7ia,+ , a € R , defined 
below. 

Definition 1.1 Ha,+ is the Hilbert space of R" x R" -valued functions of z & Z'^ endowed 
with the norm 

nii+= E \Y{z)\\i + \zrr<oc. 

In addition, it is assumed that the initial data vanish ( Fq = ) at = . 
We impose the following conditions E1-E6 on the matrix V . 

El. There are positive constants C and 7 such that ||T^(-z)|| < Ce~^^^^ for 2; e Z'^ , where 
||\^(2;)|| stands for the matrix norm. 

Let V{9) be the Fourier transform of V{z) with the convention 

zez<* 

where " ■" stands for the inner product in Euclidean space R'^ and T'' for the c? -torus 
R7(27rZ)'^. 

E2. 1/ is real and symmetric, i.e., Vik{—z) = Vki{z)&'R,, k,l = l,...,n, z & Z'^ . 

The two conditions imply that V{0) is a real-analytic Hermitian matrix-valued function 
of ^ e T<^. 

E3. The matrix V{9) is non-negative definite for every ^ e T*^ . 
Let us define the Hermitian non-negative definite matrix, 

n{e) = {v{e)f^ > o. (1.12) 

The matrix Q{9) has the eigenvalues < u;i{9) < U2{0) . . . < ijJs{9) , s < n, and the 
corresponding spectral projections with multiplicity = trUaiO) . The mapping 

9 I— > uj„{9) is the cr-th band function. There are special points in T"^ at which the bands 
cross, which means that s and jump to some other value. Away from such crossing 
points, s and are independent of 9 . More precisely, the following lemma holds. 
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Lemma 1.2 (see [7, Lemma 2.2]). Let conditions El and E2 hold. Then there exists a 
closed subset C* C T*^ such that the following assertions hold: 
(i) the Lebesgue measure of C* is zero; 

(a) for any point © e T'^\C^, , there exists a neighborhood 0{Q) such that each band function 
uJa{9) can he chosen as a real-analytic function on 0(Q) ; 
(Hi) the eigenvalue uJaid) has constant multiplicity in T'^ \ ; 
(iv) the following spectral decomposition holds: 

n{e) = j2uja{e)Ti,{e), eeT<^\c, (1.13) 

CT=1 

where ^a{(^) is an orthogonal projection in R", and is a real-analytic function on 
T'^ \ C . 

For 9 & , denote by Hess {uja-) the matrix of second partial derivatives. The next 

condition on V is as follows. 

E4. Let = det(Hess(u;^(^))) . Then Da{9) does not vanish identically on T'^ \ C* , 

Let us write 

Co^{9eT'^: det V^(^) = 0} and = {9 e \C^ : D^{9) = 0}, a = l,...,s. (1.14) 

Then the Lebesgue measure of Ca vanishes, a — 0,1, s (see [7, Lemma 2.3]). 
The last two conditions on V look as follows. 

E5. For each a ^ a' , the identities 0J(j{9') ± ijJa'{9') = const± for ^ e T*^ \ , do not hold 
with const± 7^ . 

This condition holds trivially for n — 1. 

E6. ||V^-^(^)|| e L\T'^). 

If Co = 0, then HV'-H^)!! is bounded, and E6 holds trivially. 

Remark 1.3 Conditions E1-E6 are satisfied, in particular, in the case of the nearest neigh- 
bor crystal in which the interaction matrix V{z) = {Vki{z))l^i^i is of the form 



Vkliz) =0 iOT k^l, Vkkiz) = 



-7fe for \z\ = 1, 

27fe + m^ for 2; = 0, k — l,...,n, 

for 1^1 > 2, 



with 7fe > and > . In this case, equation (1.5) becomes 

Uk{z, t) = (TfcAi - ml)uk{z, t), k^l,...,n. 
Here stands for the discrete Laplace operator on the lattice Z'' , 

Alu{z) := ^ {u{z + e) -u{z)). 

e,|e|=l 
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Therefore, the eigenvalues of V{9) are 

Uk(0) = ^J 271(1 - cos ^1) + ... + 27d(l - cos Od) + ml , k^l,...,n. (1.15) 

These eigenvalues still have to be labelled according to magnitude and degeneracy as in 
Lemma 1.2. Clearly, conditions E1-E5 hold with = . If > for any k , then the set 
Co is empty and condition E6 holds automatically. Otherwise, if nik = for some k , then 
Co = {0} . In this case, E6 is equivalent to the condition uj^'^iO) G L^(T'') , which holds if 
(i > 3 . Therefore, conditions E1-E6 hold if either (i) 0? > 3 or (ii) d = 1,2 and > for 
any k . 

Lemma 1.4 (see [10, Corollary 2.4]) Let conditions El and E2 hold. Choose some a G R . 
Then (i) for any Yq G V-a,^ , there exists a unique solution Y{t) G C(R, to the mixed 

problem (1.8); 

(ii) the operator U-^{t) : Yq i— > Y{t) is continuous on Ha,+ ■ 

The proof is based on the following formula for the solution X{t) of (1.11): 

X{t)^ E Gt{z-z')Xo{z'), (1.16) 

where the function Gt{z) has the Fourier representation 

Qtiz) := Fi_l^[exp{Aie)t)] = {^t:)-" j e"^^'^ exp (i(^)i) (1.17) 

with 

Therefore, the solution Y{t) of (1.8) admits the representation 

Y{t)^ QtA^^^'W)^ ^gZ^, (L19) 

where ^t,+ (z, z') := Qtiz - z') - Stiz - z'). (1.20) 

2 Main results 

2.1 Convergence to equilibrium 

Denote by fiQ a Borel probability measure on Ti.a,+ giving the distribution of Yq . The 

expectation with respect to ^0 is denoted by Eq . 

Assume that the initial measure fiQ has the following properties S1-S4. 

SI. ^o(^) has zero expectation value, Eo{Yq{z)) = for 2; G Z"^ . 

For a,b,c E C" , denote by a ® 6 the linear operator (a ® 6)c = a J^j=i bjCj . Denote by 
Ha the Hilbert space of R" x R" -valued functions of 2; G Z'^ endowed with the norm 

ll^lla= E \X{z)\\l+\z\'r<OC. 
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Definition 2.1 A measure v is said to he translation invariant if v{ThB) = y{B) for 
B e and heZ'^ , where ThX{z) = X{z - h) for z e Z'^ . 

52. The correlation matrices of the measure /io have the form 

Qi^(z,z')^Eo{Y^{z)®Y^(z'))^qi\z„z[,z-z'), z,z'eZ'i, i,j ^0,1, (2.1) 

where (i) {zi, z[, z) = for Zi = or z[ = , 

(ii) limy^+^qQ {zi+y,y,z) = (io (z) , z = {zi,z) e Z*^ . Here c[o (z) are correlation functions 
of some translation invariant measure uq with zero mean value on Ha ■ 

53. The measure /Iq has finite variance and finite mean energy density, 

eo{z) = Eo{\Y,\z)\' + \Yo\z)\') = tr [Q','{z,z) + Ql'{z,z)] < eo < oo, ze Z%. (2.2) 

Finally, it is assumed that the measure /xq satisfies a mixing condition. To formulate 
this condition, denote by (j(^A) , A C Z5J[_ , the a -algebra on 7^0,,+ generated by Yo{z) with 
z A. Define the Ibragimov mixing coefficient of the probability measure \xq on 7^0,+ by 
the rule (cf. [11, Definition 17.2.2]) 

/ N |/xo(^n5) -/xo(^)/^o(5)| .oQ^ 
(^(r) = sup sup — ^ — (2.3) 

dist(^, B)>r Mo(S)>0 

Definition 2.2 A measure fiQ is said to satisfy the strong uniform Ibragimov mixing condi- 
tion if (p{r) — >• as r — >• oo . 

54. The measure //q satisfies the strong uniform Ibragimov mixing condition with 

r*^- V^/^(r) dr<oo. (2.4) 



/- 



Remark 2.3 The uniform Rosenblatt mixing condition [13] is also sufficient, together with 
a higher power > 2 in the bound (2.2). Namely, there is a (5 > such that 

Eo{]Y^{z)\'+' + \Y^{z)\'^') < C < oo, zeZ%. 

+ 00 

Condition (2.4) needs a modification, namely, j r'^~^a^{r)dr < oo with p = min(5/(2 + 



5). 1/2). Here a{r) is the Rosenblatt mixing coefficient defined as in (2.3) but without 
/io(i?) in the denominator. Under these modifications, the statements of Theorem 2.5 and 
their proofs remain essentially unchanged. 

The uniform Rosenblatt mixing condition can also be weakened, see Remarks 3.4 in [10]. 

Definition 2.4 We define fit as the Borel probability measure on Ha,+ which gives the 
distribution of the random solution Y{t) , 

pi^{B) = fjLo{U+{-t)B), where B e B{na,+) and teR. 
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In [10], we prove the weak convergence of the measures /it on the space 7ia,+ with 
a < —d/2 to a hmit measure /Xoo , 

fJ't — ^ fJ'oo as i ^ oo, (2.5) 
where /^oo is an equihbrium Gaussian measure on 7ia,+ ■ This means the convergence 

^Um| f{Y)f,t{dY) = J f{Y)f,^{dY) 

for any bounded continuous functional / on 7ia,+ ■ 

Theorem 2.5 (see flOj). Let d,n >l,a< —d/2. Assume that conditions (1.9), E1-E6, 

and S1-S4 hold. Then 

(i) the convergence in (2.5) holds. 

(a) The limit measure fi^o is a Gaussian measure on 'Ha.+ ■ 

(Hi) The correlation matrices of the measures fj^t converge to a limit for i,j = 0,1, 

Qi'iz, z') = j {Y\z) Y^(z')) iXtidY) ^ Q%{z, z'), t^ oo, z, z' e Z%. (2.6) 

The correlation matrix Qoo{z,z') = {Q'^Joi^j z'))i,j=o of the limit measure has the form 
Q^{z,z')=q^{z-z')-q^{z-z')-q^{z-z') + q^{z-z'), z^z'eZ'^. (2.7) 

Here qoc,{z) = q^{z) + q^{z) , where in the Fourier transform, we have 

s 



qM - lj:^A0){U0) + cmo{0)c{9r)n,{e), (2.8) 

^ a=l 

= li2sign(do,u,(e))u,(e){c(e)cio(e)-cio(o)cm^.(o), ^ e T'^\a,(2.9) 

^ <T=1 

lla{0) is the spectral projection in Lemma 1.2 (iv), and 

(iv) The measure fi^ is time stationary, i.e., [t/+(t)]*/Xoo = A*oo ; t € R- 

Remark 2.6 (i) From formulas (2.8)-(2.10) it follows that qoo{Q) satisfies the "equilibrium 
condition", i.e., one has the form 

a (f))-( ^ (2n) 

Qoo[t^) - \^ n'^(e)h(e) ) ■ ^^-^^^ 

Moreover, {q%{e)r = q^{e) , ^,j=0,l,and > , q'M > . 

(ii) Introduce the complex-valued field 

a{x) = -^{yl^\o{x) + iV+^^\i{x)) e C" , X e (2.12) 
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where 

V> := Yl i^'^i^ -^)- ^""(^ - ^)H^)^ with V''{z) := Fg_^^ {v^{0)) ■ 



Let a{x)* stand for complex conjugate field. Obviously Et{a{x)) — 0. The covariance 
Qoo{x, y) has two parts. By Theorem 2.5, the aa-, equivalently the a*a* -, covariance satisfies 

lim EAaix) ® aiy)) = 0. 

For the a* a -covariance, Theorem 2.5 (ni) imphes 

lim EAaixY ® aiy)) = Wix - y) - Wix - y) - Wix - y) + W(x -y), x,y e Z+, 

t— »oo 

where in Fourier transform 

w{0) = m fM + iCid) = ^ E (i - sign{Wo,u;M))u,{e)Womaie), 

with Wo{e) = I MO)q^'{9)+u;-\0)ql\e)+iQ'o\O) -iq^V)] ■ 

2.2 Initial measure with slow variation 

Let {fJ.Q,e > 0} be a family of initial measures. Roughly, in a linear region of size , 
£ <S 1 , ijLq looks like the initial measure from Section 2.1. To formulate the main conditions 
V1-V2 on the initial covariance, let us introduce the complex 2nx2n matrix- valued function 
R{r, X, y) — {K^^lr, x, y))^j^Q , r e R*^ , x,y & Z'^ , with the following properties. 

10. R{r,x,y) = for xi = or yi = . The nx n matrix- valued functions R^^{r,x,y) 
have the form 

i?'^(r,x,y) = R'^(r,xi,yi,x - y), where x^{xi,x), y = ij^ 0,1. 

Moreover, 

lim R'^{r,yi + zi,yi,z)^Ri'{r,z), z ^ {z,, z) e Z'', i,j^O,l. (2.13) 

jyi— >+oo 

11. For every fixed r e R'^ and i, j = 0, 1 , the bound holds, 

\R^{r,x,y)\<C{l + \x-y\)-\ x,y ^Z%, (2.14) 
where C is some positive constant, 7 > d . In particular, for every r e R*^ , 

\^i{r,z)\<C{\^\z\)-\ zeZ''. (2.15) 

12. For every fixed r e R'^ , the matrix- valued function it! satisfies 

i?^^(r,-,-)>0, R^{r,x,y) = {R^\r,y,x))^, x,yeZ%. 
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lC^{r,e) (2.17) 



In particular, for every fixed r G R'^ , Ro('"i^) satisfies 

R°°(r, 6) > 0, Kl\r, 6) > 0, R°^(r, 6) = R^°(r, ^)*, Oe T". (2.16) 

13. For every fixed r G R*^ and ^ G T'^ , the matrix Ro(r, ^) is nonnegative definite. 

14. For every ^ G T*^ , Rq (•, ^) , i, j = 0, 1 , are fiinctions and the function 

Qa 

r sup max max — — . 

0^T.di,j=O,la={a,...,ad):\aj\<l . . . OT^ 

is bounded uniformly on bounded sets. 

Remcirk 2.7 For simphcity of proof, we could assume that R{r, x, y) has the simpler form, 
namely, 

i?(r,x,y) = C(.Ti)C(yi)Ro(r,a;-i/), (2.18) 

where C^{x) , a; G , is a nonnegative bounded function such that C,{x) = for x < and 
C{x) = 1 for X > a with some a > 1, and Ro(r, a;) satisfies conditions (2.15)~(2.17). Then 
R{r,x,y) satisfies 11-14. However, the limit covariance Qoo{x,y) in (2.7) has not the form 
(2.18), in general. Formula (2.7) implies that Qoo{x,y) satisfies the bound similar to (2.13), 



Qoo{x,y) = (^oo{xi,yi,x - y), lim (^ooivi + Zi,yi, z) = q^{z) + q^{z), zeZ'^. 

Therefore, we will prove the main results under condition (2.13) which is weaker than (2.18). 

Let Eq stand for expectation with respect to the measure //q . Assume that Eq{Y^{x)) — 
and define the covariance 

Q'i{x,x')^El{Y\x)®Y^{x')), x,x' i,j=0,l. 

Definition 2.8 We call a family of measures {/Xq, £ > 0} a family of slow variation for R 
if {Q*''(x, x'), £ > 0} satisfies conditions V1-V3 listed below. 
VI. For any e > , there exists an even integer such that 

(i) for all M G R'' and x, G /m n Z:j_ , then 

Q^{x,x')- R^{eM,x,x')\ < Cmin[(l + |x - x'|)-^, sA^^], (2.19) 

where C, 7 are the constants from (2.14), and Im is the cube centered at the point M 
with edge length , 

Im = {x= (xi, . . . , Xd) G Z'^ : \xj - M, | < N,/2, M = (Mi, . . . , M^)}. (2.20) 

(n) N,^e-P as e , with some p G (1/2, 1) . 

V2. For any £ > and aU x,x' G Z% , i,] = 0, 1 , \Q^^{x,x')\ < C(l + |a; - x'\)-^ with 
constants C, 7 as in (2.14). 

To prove the weak convergence of the measures (Theorem 2.14 below) we need the stronger 
condition V3: 

V3. The measures //q satisfy the Ibragimov mixing condition S4 (see section 2.1) with the 
mixing coefficients (f^ . Moreover, it is assumed that 

sup |<^y^(x)| < C(l + \x\)-^, with a 7 > d. 

£>0 

Note that condition V3 implies V2. 
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Definition 2.9 (i) fil is a Borel probability measure on 7ia+ which gives the distribution 
of Y{t) , 

lil{B) = nl{U+{-t)B), where B e B{na,+) and teR. 
(a) The correlation functions of the measure /if are defined by 

QU^,y)- I {Y'{x)^Y^{y))f,l{dY)^E',{Y'{x,t)^Y^{y,t)), i,j=0,l, x,yeZ'i. 

Here t) are the components of the random solution Y{t) = {Y^{-,t),Y^{-,t)) to the 

problem (1.8). 

2.3 Covariance in the kinetic scaling limit 

Let us use a time span of order r/e" , r ^ , with < « < 1 , and study the asymptotics of 
the covariance Qcr/e'^i^i u) s& £ — > +0 . For < a < 1 , the result is given by Theorem 2.15. 

To formulate the result for a — 1 let us introduce the matrix qr,r{z) , z e , r e R'^ , 
T 7^ , by the Fourier transform, qr,r{z) — -^e^z[?r,r(^)] ; where 

= Ena(^)[M;(r;r,e) +iM^(r;r,^)]n,(e), 9 e T''\C.. (2.21) 

a=l 



(2.22) 



Here Ila{d) is the spectral projection introduced in Lemma 1.2 (iv), 

M-(r;r,^) = i(R^(r; r, ^) + C(^)R^(r; r, ^)C*(^)), 
M-(r;r,^) = i(C(^)R- (r; r, ^) - R- (r; r, ^)C*(^)), 

with C{0) defined in (2.10), and 

R^(r; r, 9) = ^(Ro(r + Vu;,{9)t, 9)xt,r,{0) ± Mr - V^.(^)r, ^)x;,.,(^)), (2-23) 

where 

xtrM = ^(1 + sign(ri ± tV.u^M). (2.24) 

Theorem 2.10 Let conditions V1-V2 and E1-E6 hold. Then for any r 7^ , r G R'^ with 
ri > , the correlation functions converge to a limit, 

lim g,,,/,([r/£] + z, [r/e] + z') = Q%{z, z'), (2.25) 

where 

qG ^ ( <1?A^ - z') = qT,r{z - z') + qr^z - z'), ifri > 0, z, z' G Z'^; 

' \ qrA^-z') - qr,r{z-z') - qr,r{^-z') + qr,r{^-z'), if ri = 0, z,z' e Z^, 

with qr,riz) defined by (2.21)-(2.23). 
This theorem is proved in Section 3. 
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Corollary 2.11 Let n > and r 7^ . (i) From formulas (2.21)-(2.23) it follows that the 
a -hand of q^^ri^) satisfies the following "hydrodynamic" equation: 

drf{T,r-e)^iC,{e)Vu,{e)-Vrf{T,r-e), ri>0, T>0, 



where Cc — 




, a — 1, . . . , s, the boundary and initial conditions are given by 



the R-o • 

(a) From formulas (2.21)-(2.23) it follows that q^ ^. satisfi.es the equilibrium condition (2.11), 
i.e., q^'^\ey = n\e)q^^^\e), q^'^\e) = -q^'^\e). Moreover, gff(^) > 0, q^fiO)* = 

q?f{e). 

Let us introduce the scaled n x n Wigner matrix through 

W^{T;r,e)= Y.e''-'E^^/^(a*{[r/e + y/2])(^a{[r/e-y/2])), r G 

where a{x) is given in (2.12). By conditions VI and V2, the following limit exists 
limW^^(0;r,^) = ^{n'/''B^ir,9)Q'/^ + n-'/^Rl\r,9)n-^/^ 

= W^(0;r,^), ri>0. 
We also define the hmit Wigner matrix as follows. 

cr=l 

+W{0; -n + rViu;,{e),f - rVu;,{e), ^~)x+_,,(^)}n,(^) (2.26) 
t Il^{e)W(0;r-rVu;,{e),e)U,(e), if n > rViU^ie), 

E U,{e)W{O;^, + TS7e,u;,{e),f-TVeiVa{0),ma{O), if n<TViu;,{e), 

K a=l 

where § = {-d^,d) , 6 = {62, ...,e,),r = {r^, . . . ,r,) . 

Theorem 2.12 Let conditions W-Y 2 andEl-EG hold. Then for any r e and T>0, 
the following limit exists in the sense of distributions, 

lim W'{t; r, 9) = WP{t; r, 6) . (2.27) 
In addition, for the remaining part of the covariance, 

lim J2 e''-'E'^/^(a{[r/e + y/2])0a{[r/e-y/2]))^O. 

yeZd 

This theorem is proved in Section 4. 
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Corollary 2.13 Denote by W^{t; r,9) , a = 1, . . . , s , the a -th hand of the Wigner function 
W^ir^r^O) . Then is a solution of the "energy transport" equation 

drW^{T;r,9) + VLU,ie) -VrW^iT^r^e) =0, r > 0, r G 

where the boundary and initial conditions are given by the initial Wigner matrix projected 
onto the a -th band, 

W^(T;r,e)\,=o = U,{e)W{0;r,e)U^id), rGRt 
W^{T;r,9)\r,=o = b{T;f,e), f e R'^-\ r > 0. 



Here 

b{T;f,9) :-- 



U,{9)W{0; {-TV,u;^{9),r - tVu^{9)), 9)U^{9), if V^uM < 0, 
Ii,{9)W{Q-{TV^uj,{9),f-TVuj,{9)),9)n,{9), if V^uj,{9) > 0. 



2.4 Weak convergence of measures family 

Let us consider the random field Y at the kinetic time r/e , r 7^ , and close to the spatial 
point [r/e] e Z% . Denote by T/^ , /i e , the group of space translations. The measure at 
r/e is then defined through 

/^r/£,r = T^[r/e]lJ'%/e > (2-28) 

i.e., iJ,%i^,^{B) = ii%i^{T\r/e\B) , where B e B{Ti.a,+) and /x^/^ is defined in Definition 2.9. 

Theorem 2.14 Let conditions V1-V3 and E1-E6 hold. Then for r 7^ 0, r G R"^ with 
ri >0 , in the sense of weak convergence on 'Ha,+ , 

Jim/^^M. = /^?r- (2-29) 

The measure /i^^ is a Gaussian measure on 'Ha,+ , which is invariant under the time trans- 
lation U+{t) . /i^,^ has mean zero and covariance 



Q?fiz,z')^ J iY\z)^Y^iz'))^l%idY), 



defined by Theorem 2.10. If ri > , the covariance Q^^'^^{z,z') = q^^^!^ {z — z') is determined 
through W'^{T;r,9) as 

mQT{0)^m-'Q?f{0) = l{W^{T;r,9) + W^{T;r,-9f) (2.30) 



and 

qffi9) = -q?f{9) = -'-{W^{r;r,9) - W^{T;r, -9f) . (2.31) 



Theorem 2.14 is proved in Section 5. 

Prom Theorem 2.14 we conclude that close to r/e in space and close to r/e in time 
the random field Y^{x,t) is a stationary Gaussian field. Its distribution at fixed local time 
r is given by fi^^ while in time it evolves deter ministically according Uj^{t) . In this sense 
locally in space and time the random field is stationary with statistics determined through 
the Wigner matrix at (r, r) and the microscopic dynamics, compare with (2.30), (2.31). 

Let us use a time span of order r/^" with an a G (0, 1) . In this case, change condition VI 
(ii) as follows: ~ e^^ as £ — with some jS G (q;/2, a) . Then the following result holds. 
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Theorem 2.15 Let a e (0, 1) , conditions V1-V2 and E1-E6 hold. Then for r 7^ 0, 
(i) the correlation functions of measures if^j^a^ converge to a limit, 

limg,,^/^a([r/£] + z, [r/e] + z') = Qr{z, z'), z, z' e 

where Qr{z,z') does not depend on r and has the form 

Qr{z,z') = 



qr(z-z'), ifri>0; 
qr{z - z') - qr{z - z!) - qr{z - z') + qr{z - z'), ifri = 0; 



where in Fourier space, qr{0) = Ila-{0)Mr{0)Ila-{9), with 

1 



Mr (9) 



Ro(r,e) + C{9)Ro{r,e)C*(e) , ifr^ > 0; 



Mr{e) = -\R^{r,e) + c{e)ii^{r,e)c*{e) 

+i(C(^)Ro(r, 9) - Ro(r, e)C*{e)) sign(Viu;,(^))] , i/ ri = 0. 

(ii) The measures [i^jf^ar converge weakly on the space 'Ha,+ to a limit measure fir as 
£ — >• . Moreover, fir is a Gaussian measure on Ha,+ , which is invariant under the time 
translation U+{t) , has mean zero and covariance Qr{z,z') defined above. 

We omit the proof of Theorem 2.15 since it can be proved by using the technique of 
Theorems 2.10 and 2.14. 



3 Convergence of correlation functions 

In this section we prove Theorem 2.10. Before we outhne the strategy of the proof. At 

first, we use the cutting strategy from [7]- [9] combined with some techniques from [3], where 
Theorem 2.10 has proved for the case when d = n = 1 and in the entire space (see [3, 
Theorem 3.1]). Note that in [3] it is assumed the stronger conditions on matrix V than E3, 
E4, namely, a;(6') > , and the set 

{e e [-7r,7r] : cu" {9) ^ cu'" (9) ^ 0} 

is empty. Under these conditions, in [3] the uniform asymptotics of the Green function is 
proved, 

sup \gt(x)\ <c(i + \t\y^/^ (3.1) 

This bound plays an important role in the proof of [3] . However, if n > 1 , then uUs may 
be non-smooth because of band crossing, and if d > 1 , the set where the Hessian vanishes 
does not consist of isolated points. Therefore a strong estimate as (3.1) is unlikely to be 

valid, in general. To cope with such a situation, we split Qt{x) into two summands: Gt{x) = 
Qt{x) + Qt{x) ! where Qt{x) has a support in the neighborhood of a "critical set" C C T'^ , 
and Qt{x) vanishes in the neighborhood of C . The set C includes all points ^ G T"^ either 
with a degenerate Hessian of uJaiO) , or with non-smooth uj„{9) (see formula (3.4)). We show 
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that the contribution of Qt {^) is neghgible uniformly in t (see (3.9)). Hence, it allows us to 
represent correlations functions Qe,T/e in the form: Qe,T/e = Ql t/s + r/e ' ^^^^ ^^^^ 

Qir/eM= E sU^,^')Qe{x',y')GUy,yr- 

For the remainder Ql^r/e ~ Qs,T/e ~ Qlr/e prove that Ql.j./^{x,y) = o(l) uniformly 
in r 7^ , e > and x,y G Z'^ . The last fact follows from two key observations: (i) 
mcsC = (Lemma 1.2) and (ii) the correlation quadratic form is continuous in £^ , sec 
Corollary 3.3. Up to this point we apply the "cutting strategy" from [7]-[9]. Finally, we prove 
that T-/e([^/£] +2;, [r/s] +y) converges to a limit as £ — > , using the techniques of [3] and 
[8]. In addition, the asymptotics of (x) , (see Lemma 3.4) of the form (x) ~ (1 + |t|)~''/^ 
plays the important role, since it replaces the asymptotics (3.1) and also simplifies some steps 
of the proof of [3]. However, in our case the structure of R{r,x,y) is more complex than 
in [3] or [9], in which R{r,x,y) = Ro(r, x — y) . We apply the approach of [8, 10], where 
convergence to equilibrium was proved for non translation-invariant initial measures, and 
combine with the technique of [9]. 

3.1 Bounds for initial covariance 

Definition 3.1 By & = l^i^TJ^) ® R" {hy l\ = e'{Z%) ® R") , where p > I and n > 
1, denote the space of sequences f{z) — {fi{z),...,fn{z)) endowed with norm \\f\\ip — 

( E \f{z)\^y^' , respectively, ||/||,. := ( E |/(^)h'^'- 

The following lemma follows from condition V2. 

Lemma 3.2 Let condition V2 hold. Then, for i, j = 0, 1 , the following bounds hold: 

\Q'Jiz,z')\ < C<oo forallzeZ%, 

z'ez'\_ 

Y.\Qe{^^^')\ < C<oo forallz'eZ%. 

Here the constant C does not depend on z,z' ^T/^^ and £ > . 
Corollary 3.3 By the Shur lemma, it follows from Lemma 3.2 that 

\{Q,{z,z')M^)®-^{^'))A<C\\nt\\\ne\. for any $,*e£^, 
where the constant C does not depend on £ > . 

3.2 Stationary phase method 

By (1.17) and (1.18) we see that Qt{9) is of the form 

Afa\ ( cosVtt sinOtf]-M , . 

^'^^>-\-sinQtn cosQt J' ^^-^^ 
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where Q = fl{6) is the Hermitian matrix defined by (1-12). Hence, we can rewrite Gt{x) in 
the form 

St{x)= J e-'^-^e^'^^^'^'a^{e)de, (3.3) 

by (1.13). We are going to apply the stationary phase arguments to the integral (3.3) which 
require a smoothness in 9 . Then we have to choose certain smooth branches of the functions 
a^{0) and 00^(0) and cut off all singularities. First, introduce the critical set as 

C = CoUC* U (c.U {9eT'\C.: ^^ = 0}[j{9eT''\C.: Via;,(^) = 0}^ , (3.4) 

with as in Lemma 1.2 and sets Cq and defined by (1.14). Obviously, mesC = (see 
[8, lemma 7.3]). Secondly, fix an S > and choose a finite partition of unity, 

f{0)+g{0)^l, g{0)^j29k{O), 0eT^, (3.5) 

k=l 

where f,gk are nonnegative functions in C^(T'^) , and 

supp/ C e T*^ : dist(e,C) < 5}, suppgk C {6 e T"^ : dist(^,C) > 5/2}. (3.6) 
Then we represent Qt{x) in the form Qt{x) = Qt{x) + Qfi^) ) where 

g({x) = (27:)-' J e-'^-'f (6) gt{9)de, (3.7) 
gf{x) = {27:)-' j e-'--'g{e)gt{e)de= j^j^j gk{e)e-'^-'^'^''^'^'^ (3.8) 

Td ±,0-=l fe=lrj,d 

By Lemma 1.2 and the compactness arguments, we can choose the supports of gk so small 
that the eigenvalues uJ^jiO) and the amplitudes a^iO) are real-analytic functions inside the 
supp gk for every k . (We do not label the functions by the index k to simplify the notation.) 
For the function ^/(x) , the Parseval identity, (3.2), and condition E6 imply 

WGim^-c J\g,{e)\Vm'de<c J igm' de ^ o as 5^0, (3.9) 

Td- dist{e,C)<S 

uniformly in t G R. For the function gf{x) the following lemma holds. 

Lemma 3.4 (see [9, Lemma 4-5]) Let conditions E1-E4 and E6 hold. Then the following 
bounds hold. 

(i) sup,,z^|^f(x)|<cr^/2. 

(a) For any p > 0, there exist Cp,'jg > such that \gf{x)\ < Cp{\t\ + \x\ + 1)~^ for 
\x\ > -fgt . 
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3.3 Proof of Theorem 2.10 

The representation (1.19) yields 

Q,,t{z,z')^E^,{Y{z,t)^Y{z',t))^ ^ gt,4z,x)Q,{x,y)gt,+{z',yf, z, z' e Z^, (3.10) 

for any t G . It follows from condition (1.9) and from formulas (1.17) and (1-18) that 
Gt{z) = Qt{,z) with z = {—zi,Z2, . . . ,Zd) ■ In this case, by (1.20), the covariance Q^^t{z,z') 
can be decomposed into the sum of fourth terms, 

Q,,tiz, z') = S,,tiz, z') - S,,tiz, z') - S,,t{z, z') + S,,t{z, ~z'), z, z' e 

where 

S,,t{z,z'):= Yl Qt{z-x)Q,{x,y)gt{z' -yf. (3.11) 

Proposition 3.5 Let r E R'^ and z, z' E Z'^ . Then 

Se,r/e{ir/e] + z, [r/s] + z!) ^ g,,,(^ -z!), +0, (3.12) 

where qT,r{z) is defined in (2.21)-(2.23). 

This proposition implies Theorem 2.10. Indeed, let ri — . Then f/e — r/e — (0, f/s) 
and for z, e , 

Qs,r/e([r/£] + z, [r/e] + z') = S,,r/e{[r/e] + z, [r/e] + z') - S,,r/e{[r/e] + z, [r/e] + z') 

-S,^r/e{[r/e] + z, [r/e] + z') + S,,-,,,{[r /e] + z, [r/e] + z'). 

Therefore, convergence (3.12) imphes (2.25). 

Let ri > 0. In this case, the matrix- valued functions Se,T/e{[f/^] + z, [r/e] + z') and 
Se^r/eiV/^] + z, [f/e] + z') vauish as £ ^ +0 , and 

Se,r/e{[f/e] + [f/e] + z') ^ qr,f{z -z'), £ ^ +0. 

It can be proved by similar way as Proposition 3.5. 
Proof of Proposion 3.5. Step (i) Let us denote 



Qe{x,y) for x,y E Z\, 
otherwise 



Qeix,y) = 

Corollary 3.3 and (3.9) imply that 

S,,t{z,z')^ Y: gf{z-x)Q,{x,y)g!{z'-yf + o{l), 

where o(l)— >^0as5— >0 uniformly in i e R and z,z' E Z''' . In particular, setting t — r/e , 
z — [r/e] + I and z' — [r/e] +p we get 

S.,r/.{[r/e]+lAr/e]+p) = ^ g'r/,{[r/e] + I - x)Q,{x,y)g^^/^{p+ [r/e] - yf + o{l) 

x,yeZ'^ 

= E G'r/eii + ^)Qei[r/e]-x,[r/e]-y)g^^/^{p + yr + o{l). 
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Let c = 7g + max(|/|, |p|) . Then Lemma 3.4 (ii) and condition V2 imply that 

5,,./,([r/£] +1, [r/e] + p) = ^ ^r/ei^ +x)Q.{[r/e]-x, [r/e]-y)g'^/^ip + yf 

+ri(£,r) + o(l), £^0, 
where hm£~Pri(£,T) = for any p > and r e . 

£—♦0 

Step (ii) We divide the cube [—cr/e^cT/e]''' onto the cubes InN^ (see (2.20)), 

[-CT/e,CT/e]'^ C U /nJV,, 

where J = {n = {rii, . . . .Ud) e 7/ , \nj\ < [cr/{eNs)] + 1} . Therefore, 
S,,r/,{[r/e]+l, [r/e] + p) = E E ^'^(^ +x)Qe{[r/e] -x, [r/e] -y)QliXv + vf + o(l)- 

Now we prove that the contribution of the sums over pairs m, n e J with m ^ n vanishes 
as £ — > . Let us denote 

r2{e,r) = ^ Ql,e{l +^)Qe{[r / e]-xAr / e]-y){gi,^{p + y)^ (3.13) 

m,n & J, m ^ n 

X e ImNe , y e InNs 

and prove that 

r2(£,T)^0 as £^0 (3.14) 

for any r G . We divide the sum in the RHS of (3.13) onto two sums and 5*2 , 

where the first sum 5*1 is taken over all x G ImN^ , V G InN^ and m,n & J such that 
3j G {1, . . . ,d} : \mj — nj\ > 2; the sum 5*2 is taken over all x G ImN^ , y G InN^ and 
m,n E J such that m ^ n and Vj = 1, . . . , o? : |mj — n^l < 1 . The number of points m E J 
is order of (r/ (eiV^))'^ , the number of points a; G ImNe is ~ ■ Therefore, by Lemma 3.4 
(i) and condition V2, the sum is estimated by 

C{l + r/er\r/er E (1 + kl)^^ 

seZd;\s\>Ne 

which vanishes as £ — >• , since iVg — >• +oo as £ — >• and j > d. To estimate the 
second sum S2 (the contribution of nearest neighbors ImNs a-nd InNs ) we choose a number 
p > d + 1 and divide the sum onto two sums: S2 — S21 + S22 , where the sum 5'2i is taken 

over all m G J and x G ImN^ , n E {n E J : n ^ m,\/j : \mj — %| < 1} and y G InN^ 
such that |a; — > Nl^'^ and the second sum 6*22 is taken, respectively, over y such that 
\x — y\ < N^^P . The contribution of "non-boundary zones" S'21 is 

C(l + r/£)-'^(r/£)'^ Y: {l + \s\)-^ 

seZ'i,\s\>N^^^ 

which vanishes as £ — > . The contribution of "boundary zones" 5*22 is order of 

C(l + T/e)-'^{T/eN,YNyP+'^-^N^/P ~ CiVi'^+^^/f-^ (3.15) 
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Indeed, the number of points m € J is order of {T/{eN^)Y , the number of points {n : 
l^i ~ ^jl < l,Vj, m n} is finite. For fixed m,n the number of points x e ImN^ such 
that \x — y\ < N^P is order of N^'^NI^'p . For fixed x the number of points y such that 
— y\ < A^e^^ is ~ N^/P . The number p is chosen such that {d+ l)/p — 1 < 0. Hence, 
(3.15) vanishes as £ —> by condition VI (ii). The decay (3.14) is proved. 
Therefore, 

S,,r/,{[r/e]+l, [r/e] +p) = E E ^r/si^ +x)Qe{[r/e]-x, [r/e]-y)g'^/^{p + yf + o{l). 

Step (in) Now we can apply condition VI (i) at the points [r/e] — x,[r/e] — y of the same 
cube I[r/e]-m.N^ and obtain 

\Q,{[r/e]-x, [r/e]-y) - R{e[r/e] - emN„ [r/e]-x, [r/e]-y)\ < Cmin[(l + \x- y\)-\eN,], 

where, by definition, the function R is equal to 

fs/ \ f R(r,x,y) iix.yeZ'i, 
Rir.x.y) ^ { ^ \ ' ^1 . 
^ ' '^^ \ otherwise 

Therefore, 

S,,r/e{[r/e\+lArle]+p)^ E E Ql/^l ■ ^Q'r/eiP + vf + o{ll (3.16) 

where R{. . .) = R{e[r/e] — smN^, [r/e] — x, [r/e] — y) . Indeed, for fixed x e ImNe i 
E min[(l + |x-y|)-T,£7V,]= + E {l + \x-y\)-'^ 

Y: eN,+ E {l + lsir^^ieN^y-'/^ 

5:(l+|s|)<(Af,)-l/7 5:(l + |s|)>(iVe)-V7 



By Lemma 3.4 (i) and condition VI (ii), we obtain 



E E Keii+^)\^^[ii + \^-y\r,^miG'r/Ap + y)r\ 

me J x,yeImNs 

< C(l + T/e)-\r/{eN,))'^N^{eN,y-'^/'' ~ ^{i-m-d/i) ^ q, e ^ 0, 



since P < 1 and 7 > d . 

By the similar arguments as in steps (i) and (ii) of the proof, the sums in the RHS of 
(3.16) can be taken over {m e J, x e ImNe, U £ Z''} ) i-e., 

5,,./,([r/6]+Ur/£]+p)=E E Y.Ql/eil+^)Ri---)Ql,eiP + yf + o{l). 

meJ xeimNg j/ez<^ 

Step (iv) Let us split the function R into the following three matrix functions: 

R^{r,x,y) := ^Iio{r,x-y), (3.17) 

R~{r,x,y) := ^Ro(r, x - y) sign(yi), (3.18) 
R^{r,x,y) := R{r,x,y) — R^{r,x,y) — R^{r,x,y). (3.19) 
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Next, introduce the matrices 

Slr/e^Sl^/ei\rle]+lArle]+p)^Y. E E ^.'/e^ +^)^"(- ■ O^.^/elP + (3-20) 

for each a — {+, — , 0} and spht S^^r/e into three terms, 3^^^/^ — S^^^^ + S~^^^ + S^^^^ . The 
convergence (3.12) results now from the following three lemmas. 

Lemma 3.6 limS^^^^{[r /e]+l, [r/e]+p) = |-s+^(/ — p) , l,P G Z"' , where is defined as 

in (2.21)-(2.22) hut with ^(Ro(r + Vu^{e)T, 9) ± Ro(r - Vu;^(^)t, 9)) instead of R^ (cf 
(2.23)). 

Proof. By (3.17) and (3.20), the function 5"^^^^ can be represented as 

Str/e = ^ E E Q'r/eil + ^) E ^o{e[r / s] - smN,, y - x)QIiXp + yf. 
Using Fourier transform and the Parseval equality we can rewrite S^^^^ as 

Str/e = (2vr)"^^iE E /e-('-^---^\--(^'-^)^,^/,(^)Ro(e[r/.]-.miV,,^') 

xQll,{9'yd9'd9. (3.21) 

Therefore, the proof of Lemma 3.6 reduces to the finding the limit value of (3.21), that is 
done in Theorem 4.1 from [9] (the detailed proof see in Appendix A). 

Lemma 3.7 liniS'~^y^([r /£] + /, [r/e\-\-p) — \s~,^{l — p) , l,p G Z'^, where s'^. is defined as 
in (2.21)-(2.22j hut with 

^(Ro(r + Vuj,{9)t, 9) sign(ri + tV^uj,{9)) ± Ro(r - Vu;,(^)t, 9) sign(ri - tV^uj,{9))) 
instead of R^ . 

Lemma 3.8 lim5°^/^([r/£] + /, [r/e\+p) = 0, l,peZ'^. 

The proofs of Lemmas 3.7 and 3.8 see in Appendices B and C, resp. 



4 Convergence of Wigner matrices 

Here we prove Theorem 2.12. Theorem 2.10 implies that for any r £ R+ , r 7^ and 
y e (2Z)'^ , the following convergence holds, 

hmE^^/M[r/e]+y/2y a{[r/e] - y/2)) = >V^(r;r,y), (4.1) 
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where in the Fourier space one has 

= WP{T;r,e), (4.2) 

by formulas (2.26) and (2.21)-(2.23). Therefore, convergence (2.27) follows from (4.1), (4.2) 
and Lemma 4.1. 

Lemma 4.1 Let conditions V2 and El-ES, E6 hold and let a < —d/2 . Then the following 
bound holds: 

sup sup ||Q£,t(;2;, ;2;')ll < < oo. 
teR z,z'ez| 

Proof The representation (1.19) gives 

QU^,z') = E^,(Y^{z,t)^Y^{z',t))^ J: Y: Ol^{z,y)Q'J{y,y')gi;4^',y') 

y,y'£Z'\_ k,l=0,l 

= {Q,{y,y'),^l{y,t)®^i,{y\t))+, 
where ^l{y,t) is given by 

Kiy,t) = (gf;4^,y),gil4^,yy 



= {gr{z-y)-gf{z-y),gi\z-y)-gi\z-y)), t = 0,l. 

The Parseval identity, formula (3.2), and condition E6 imply that 

\m;mi = i^^r' I me,t)\ue<c j {\gfm' + \Qfm') de<c,<^. 

Corollary 3.3 gives now 

mt{z.z')\ = m{y.y').^\{y.t) ® ^Ay'.t))+\ < cim-Mii WA-Mii <c^<oc, 

where the constant Ci does not depend on 2;' e , i e R , and £ > . ■ 

5 Proof of Theorem 2.14 

Theorem 2.14 follows from Propositions 5.1 and 5.2. Proposition 5.1 ensures the existence 
of the limit measures of the family {/i^^^^, £ > 0} , while Proposition 5.2 provides the 
uniqueness. 

Proposition 5.1 Let conditions V2 and E1-E3, E6 hold. Then for any r G R'^ with 
> , r 7^ 0, the family of measures {iJ^'^/^^j.^ e > 0} is weakly compact on 7ia,+ for any 
a < —d/2 , and the following bounds hold: 



sup / \\Yo\\li^l/,^,{dyo) < oo. (5.1) 
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Proof. Definition 1.1 yields 

/ ||lo||^,+/^^/.,.W = £^o^(||T-[./.]C/+(r/£)yo||^,+) 
= E (1 + \A'nAQZ/e{[r/e\ + [r/e] + z) + g^,V/,([r/e] + [r/e] + ;.)). (5.2) 

Since a < —d/2 , estimate (5.1) follows from Lemma 4.1 and (5.2). Now the compactness of 
the measures family {/x^/g ^ > 0} follows from the bound (5.1) by the Prokhorov compact- 
ness theorem [16, Lemma n.3.1] by using a method applied in [16, Theorem XII. 5. 2] because 
the embedding l~La,-\- C 'Hj3,+ is compact for a > /3 . ■ 

Set = [S'(Z^) X R"]^ . where 5'(Z':|_) stands for the space of rapidly decreasing real 
sequences on Z% . Write = + for Y = {Y",Y^) e 7^^,+ and 

^ = (^0, ^1) e S+ , where {Y\ = Ezez^ Y'i^) ' *'(^) , i = 0,l. 

Proposition 5.2 Let conditions V1-V4 and E1-E6 hold. Then for any r e R*^ with 
ri > , T 7^ and ^ e <S+ , the characteristic Junctionals converge to a Gaussian one, 

-■= j e'^'^'^'^^l^r/eAdY) ^ exp{ - ^Q.,.(*,^)} =: f^^^ as e ^ 0, (5.3) 
where Qr,r is the quadratic form with the matrix kernel iQr,ri^Ty))i,j=o,i , 

Q.,,(*,*)= Y: E {Q?fiz,z'),^\z)®<b\z')), *ecS+. (5.4) 

i,j=0,l z,z'&Z'l 

To prove Theorem 2.14 it remains to check Proposition 5.2. Let us rewrite (5.3) as 

/i^/,,,(*) =£;o^(exp{i(r_[,/,]C/+(T/£)yo,*)+}) ^/i?.(*), e^O. (5.5) 

We derive (5.5) by using the explicit representation (1.19) of the solution Y{t) , the Bernstein 
'room-corridor' technique and the approach of [9, 10]. The approach gives a representation of 
{T_[r/e]U+{T/e)YQ,'if)^ as a sum of weakly dependent random variables (see formula (5.12) 
below). In this case, (5.5) follows from the central limit theorem under a Lindeberg-type 
condition. 

5.1 Duality arguments 

In this section, we evaluate the inner product {T-[r/£]U-^-{T/e)YQ, . Introduce the function 
^*(^) as ^*(2;) = 'if{z) for 2:1 > , and 'if^{z) = otherwise. Therefore, 

{T_[r/e]U+{T/e)Yo,^z))+ = {Yo{z'),M^',T/e))+, (5.6) 

where, by definition, the function ^r{z',r/e) is equal to 

Mz',r/e) = E Gr/eA^,z'mr/e]'^4z) 

= (27r)-'' J{e-'''-^ - e-'''-^)g;{e)e'^'/'^-^^,{e) dO. (5.7) 
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Let us denote 

5° = e 5 = [^(Z'') ® R"]2 : ^{e) = in a neighborhood of C}, (5.8) 

where C is defined in (3.4). Since mes C = it suffices to prove (5.5) for ^'^ e 5" only. For 
the function , the foUowing lemma holds. 

Lemma 5.3 (cf. Lemma 6.3 from [10], Lemma 5.2 from [9]). Let conditions E1-E4 and 
E6 hold. Then, for any chosen ^* e 5° , the following bounds hold. 

(i) sup,gz^|$,(2;,r/£)| <C£'^/^ 

(ii) For any p > , there exist Cp > and 7 = 7(r, r) > such that 

\^r{z,T/e)\ < Cp{l + \z\+T/e)-P, \z\ > -fT/e. (5.9) 
This lemma follows from (5.7), (5.8), (3.2), and the standard stationary phase method. 



5.2 Bernstein's room-corridor' partition 

Write t = r/e . Let us introduce a 'room-corridor' partition of the half-ball {z G Z^j_ : \z\ < 
jt} with 7 in (5.9). For t > , we choose At and pt e N . Choose a 5 , < S < 1 , and 

Write /it = At + pt and = j/tt , 6^' = a^' + Aj , j = 0, 1, 2, . . . , nt = [(7^)//^*] . We refer 
to the slabs ^ {z e Z'^ : \z\ < Ufht, < Zi < V} as the 'rooms', to C/ = e : 
\z\ < Ufht, V < Zi < a^^^} as the 'corridors', and to Lf — {z Z^ : \z\ > Utht} as the 'tail'. 
Here z = {zi, . . . , z^) , At stands for the width of the room, and pt for that of the corridor. 
Denote by Xt the indicator of the room Rl , by that of the corridor , and by rjt that 
of the tail Lt . In this case, 

E,[xK^) + e/(^)] + ^t(^) = i, zezi, 

nt — 1 

where the symbol J2i stands for the sum ■ Hence, we obtain the following Bernstein-type 

i=o 

representation: 

Introduce the random variables rl , cl , If hy 

ri = (Fo, X^$r(-, t))+, 4 = (Fo, e/$.(-, t))+, k = (Fo, ?7t$.(-, t))+. 
Therefore, it follows from (5.6) and (5.11) that 

{T_^r/e]U+{t)Y,, = (Fo, $.(•, i))+ = E,(^l + c^') + k. (5.12) 



22 



Lemma 5.4 Let conditions VI and V2 hold and e 5° . The following bounds hold for 
t>l: 



E'o\ri\' < C{^) A,/t, Vj, 
E'oWtl' < cm pt/t, Vj, 

?e\l |2 



(5.13) 

(5.14) 
(5.15) 



E'o\k\' < Cp{^)t-^, Vp>0. 

The proof is based on Lemmas 3.2 and 5.3 (see also [9, Lemma 7.1]). 

Further, to prove (5.5), we use a version of the central limit theorem developed by Ibrag- 
imov and Linnik. If Q^_r(^, ^) = , then the convergence (5.3) is obvious. Indeed, then, 

Ea( cxp{z(yo, $.(■, r/s))+}) - A^,(M/)| = E^il cxpO(yo, r/e))+} - 1|) 

1/2 



= {{Q,{x, y), $,(x, t/e) ® ^r{y, r/e))+f' = (Q.,r/e,r(*, *))'/', 



(5.16) 



where Qe,T/£,r(*, *) ^ QT,r(*, *) = , £ ^ . Therefore, (5.3) follows from Theorem 2.10. 
Thus, we may assume that, for a given e <S° , 



Lemma 5.5 The following limit holds, 



nt 



(f)''%(i+Ar]+<-o, t^oc. 



2pt 



t 



(5.17) 

(5.18) 
logi. Then 



Indeed, (5.10) imphes that ht — pt + \ '-^ 7-— , i — > 00 . Therefore, ~ ~ 

log t ht 

(5.18) follows by (5.10). 

For simplicity, we put t — t/s . By the triangle inequality, 

^o^(exp{^(yo,$.(-,t))+}) -/i^.(^)| < |i?o'(expO(yo,$r(-,t))+}) -Eo'(exp{z5:TO) 

e^P{ - ^E,^o(IWI')}-exp{ - k.,.(*,*)} 



+ 

+ 



i^o^(exp{^E/a)-exp{ - ^E.^od^*!')} 
/i + /2 + /a. 



(5.19) 



Wc arc going to show that all summands /i , I21 I3 tend to zero as i — > 00 
(i) Eqn (5.12) implies 



7i = |£;o^(exp{i$:.rO(exp{i$:.4 + i;,}-l)) 

< E/o(l^l) + ^o(l^*l) < E,(^o(l^l'))'^' + (^1 

From (5.20), (5.14), (5.15) and (5.18) we obtain that 

h < Cnt^pt/tf^ + Cpt-P ^0, t^oo. 



1/2 



(5.20) 
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Step (a) By the triangle inequality, 



1 



hi 



+ 



(5.21) 



where Qe,t,r is the quadratic form with the matrix kernel Q^e,t,r{-'^i ^) ■ Theorem 2.12 implies 
that l2i\t=rje ^ £ — > . As for 722 , we first obtain that 



/22 < E \El{rlr\) 

3<l 



(5.22) 



The distance between the different rooms Rl is greater or equal to pt . Then, by Lemma 5.3 
(i) and condition V2, 

3<l 

< Ct-'j:i:xi{x)j:xl{y){i + \x-y\)-^ 

3<l X y 

+00 

~ t-'^ny-^At J (1 + s)-^s'^-^ ds - nt{l + Pt)"^+^ (5.23) 

pt 

which vanishes as t ^ 00 because of (5.18) and j > d. Finally, it remains to check that 
I23 , t —>■ 00 . We have 

Q,,,.(*,*) = £;o^((yo,$.(-,i))^) = £;o^((E,(^t + ^) + ^*)')' 

according to (5.12). Therefore, by the Cauchy-Schwarz inequality, 

I23 < |^o^((E/t)')-^o((E/* + E,^ + y')| 

< 



Cn,j:^E^,{\4\') + C^i(i?o^((E,W)'))'^'(^tE/o(Nf) + ^0^" 



1/2 



(5.24) 



Then (5.13), (5.22) and (5.23) imply 

mjirir) < E,^o(iw'n+2E,<J^o(^M) 

< CntAt/t + Cmtil + pt)-^+'' < C2 < 00. 
Now (5.14), (5.15), (5.24), and (5.18) yield 

I23 < Cm^tPt/t + C2nt{pt/tY''' + C^t-P ^0, t^oo. 
So, the terms I21 , /22 , I23 in (5.21) tend to zero. Then (5.21) implies that for t = t/s 

1 



l2< 



E.^o(k^r)-Qr.^(*,*) 



0, £ ^ 0. 



(5.25) 
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Step (in) It remains to verify that for t = t/s 



^o(expOE/0) -exp{ - \Y.,mri\% - 0, £ ^ 0. 



Condition V3 yields 



nt-1 



< 



E^,{eM^ri})- n i^o'(expK}) 

-nt 

nt-1 nt-1 

Eo^(exp{^r,-"'}exp{z ^ r^) - ^o'( exp{zr-"'})Eo^( exp{z ^ rl}) 

-nt+1 -nt+1 

nt-1 "t-l 



+ Eo^(exp{zrr'})^o'(exp{z ^ r^) - n^o(expK}) 

—nt+1 —nt 

nt — 1 nt — 1 

< C{l + pt)-'^+\E^,{exp{t Y: ri})- [J ^o(expK}) 

-nt+1 -nt+1 

We then apply condition V3 recursively and obtain, according to Lemma 5.5, 



nt-1 



^o(exp{^E/i})- n £^o(expK}) 



t=T/e 



<Cnr/e{l + Pr/e)-'"'^0, S ^ 0. 



Hence, it remains to show that for t — t/s 

nt—1 



n ^o^exp {iri} - exp [-^Y.,ElWt?] 



0, £^0. 



According to the standard statement of the central limit theorem (see, e.g. [12, Theorem 
4.7]), it suffices to verify the Lindeberg condition: 



1 



V5>0, -Y^^El 



CTt 



t=T/e 



as £ — > 0. 



Here at = Z^j-E'olWI > ^^'^ -^o'"/ = -^o("^"/) ' where X" is the indicator of the event 
I/I > . Note that (2.6) and (5.17) imply that ar/e ^ Qr,r(*, *) as £ ^ 0. Hence, it 
remains to verify the limit relation 

J2.Eo'yi/f ^ as £ ^ for any a > 0. 
This condition can be proved by using the technique of [9]. ■ 



6 Appendix A: Proof of Lemma 3.6 

Step (i). Let us study the sum in (3.21) over ImN^ = {x & 7/^ : {rrij — 1/2)N^ < Xj < 



(mj + l/2)N„j^l,...d}: 



ix-{0'-0) 



n 



W.-^,-,iV„m,) 



\ ^ie'r<^^ _ 1 
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where F{6j, Ns,mj) := exp{i9jNi;{mj + 1/2)} — exp{i9jNs{'mj — 1/2)} . Define the function 
iz 

a{z) : a(z) = — if 2; G (— tt, tt) \ and a(0) = 1. Changing variables in (3.21): 

{e, 9') -^{9-z,9), we obtain 



x6;^/^(^ - z)Ro{^r,e,m,e)g',/M*d^d^^ (6-1) 

where Kr,£,TO e[r/e] —smN^ . Let C{9) be defined by (2.10) and / be the identity matrix. 
Then 

g9[9) = g{9)(^cosn{9)t I + sin n{9)tC (9)), 
by (3.2) and (3.8). Let us define 

gU(^) = 9{9){cosu;,(9)tI + smuj,{9)tC^{9)) with C^{9) = (| "^^^^^^ j . 

Hence, applying the projections Ila{9) from Lemma 1.2, we rewrite the product of matrices 
in the integrand from (6.1) as (for t — r/e) 



g!{9 - z)Ko{Kr,e,m, 0)0! {9)* = E ".(^ - z)gl,{9 - z)Ko{Kr,e,m, ^)^,V l^)*".' (^) 

a,a'=l 

= ± n.(^-.)^(^-.)(Ee±^-(^-^)*^^^%^^) 
xRo(/.,,,,„,^)^(^)(Ee^^^'"^'^*^^^^^n.'W)- (6.2) 



Step (a). Let us consider the one of the terms in (6.1) (denote it by If ). The proof for the 
remaining terms is similar. 

TT 

if := (27r)"^'^- E / ^-ia-p)-e^»a>^,(e)r/e^^g^)/ f ^izj M^d)F{zd, N^, ^d) ^ ^ ^ 



^, "I ^iz2i2 ^i^2)F{z2,N„m2) / J a{zi)F{zi, N„ mQ 



— TT — TT 

xg{9 - z)U^{9 - z)Ko{Kr,e,m, 9)U^i9)dzi)dz2) . . . dzd)d9. (6.3) 

Let us write z/i = z/i(^i, ^2 - -22, • • •) = -^S^^^i^i, ^2 - -22, • • 0^/ (^^s)] > = ^^2(6*1, ^2, ^3 - 
^3, . . .) = ±[V2Cc'a(^^i, 92, 9s - Z3, ...)^/(^Ar,)] Ud = M9) = ±[V,cu.(^)t/(£7V,)] . The first 
step in the evaluating the hmit value of if is the following assertion. 

Proposition 6.1 Let condition 14 hold. Then 

TV 

;1 f ^-ie-(l-v)Au^.,ie)T/e^/n\f V- f AzM^i^'i) ^ i^'J- ■ '" '<) 



If = (27r)-2'^^ J e-'^<^-P^e'^^'^^^^''g{9)[ E / 
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\ra2-V2\<2_^ |m,i-j^i|<2_^ 
Xg±ia;.(e-.)r/en^(^ _ z)B^^Hr^^^^,e)Il,,{e)g{e - z)dz^)dZ2) . . . dZd)de + Or{l), (6.4) 

where o^(l) ^0 as £ ^ /or an?/ r 7^ . 

Proof. We generalize the strategy of the proof of Proposition 3.6 from [3], where this 
assertion is proved for d — 1 . To prove the asymptotics (6.4), we will show that the series 
in (6.3) over max^ \mj — i^jl > 3 vanishes as £ — > . 

Write Jo = {n G : \n\ < [cT/{eNe)] + 1} and hj = rrij — uj . Note that in integrand 
from (6.3) the elements of matrix product have of the form 

d 

(lla{9-z)'Ro{Kr,e,m,0)'n.a'{0)^^ = ^'T,<yl3i^ - ^)^0,l3^{l^r,e,m, 0)lla' ,^s{0) , «, 7 = 1, . . . , 0?. 

a, /3, 7, 5=1 

For simplicity of exposition, we omit the sum over a, (3, 7, 6 and assume that d = 2 . Let us 
denote M^) = 9{0)n^,ap{0) , M^) = g{9)Il^>,^s{0) , and 7^(/^l, /i2, ^) = Ro,p^{e[ri/e]- {iyi + 
hijsN^, e[r2/e] — (1/2 + h2)eNs, 9) . Hence, instead of If we evaluate the following integral: 



IZ2 



X ( j e"l(^l+^'-O ^"(^^)^C^■-^'^■^'l) e±^^^(g-^W^0,(g-;,)7^(/,l, /i^, 9)dz,)dz2)d9 . 



Here we use the fact that F{zj, N^, rrij) = e'^^^'^^'^^ F{zj, A^^, hj) . Decompose the series over hi 
and /i2 in into the sums: over hj < —3 , over \hj\ < 2 and hj > 3 , j = 1,2 . Therefore, 

= E (6.5) 

where 

/^^= E E---,4^'^^- E E---,4^'^^= E E 

/ii e Jo — 1^1 /i2 s Jo — /ii e Jo — j^i /i2 s Jo — !^2 /ii e Jo — I'l /j2 e Jo — i^2 

fti < -3 h2 < -3 /ii < -3 |/i2| < 2 fti < -3 /i2 > 3 

and so on. We want to prove that the series in 7^ over max^ \ hj\ > 3 vanishes as £ — > , 
i.e., I^'^'^^ vanish as £ — > if i = 1, 3 or j = 1, 3 . We prove this fact only for I^^'^^ . For 
remaining integrals the proof is similar. Let us write 

^min — niinj/ij e : hj < —3 and hj e Jq — Uj}. 

Hence, 



. </l2<-3_7r 
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where 



h = Il{9,Z2,h2) 

TT 

- s / 

hL,<hi<-3-^ 



izi 



(t)i{e-z) dzin{hi,h2,9).{6.6) 



To rewrite the sums over hi and /i2 we use the following "discrete integration-by-parts 
formula" (see [3, p.594]) 



E 



^i{h+l/2)NeZ _ gi(ft-l/2)iV,2 



fih) 

i(h+l/2)NeZ _ i{h 

min 



+ E [e---— -e 



/(-3) 

{f{h)-f{h + l)). (6.7) 



Let us apply (6.7) to the sum over /i2 : 

^ F{z2,N„h2)Il{e,Z2,h2) 
hl,,„<h,<-3 

1,2 



^-i5/2N,Z2 _ Qi{hl,i„-l/2)NsZ2 



Il{0,Z2,-3) 



{liie,Z2,h2)-Ii{9,Z2,h2 + l)). (6.8) 



. <h2<~4: 

mm— — 



Hence, 

|/^'^|<Ci sup sup {\Ii{e,Z2,-3)\+ \h{0,Z2,h2)-Ii{e,Z2,h2 + l)\).{6.9) 

6le[-7r,7r]2z2e[-7r,7r] . <h2<-4 

Applying the formula (6.7) to the sum over hi , we rewrite Ii in the form 

h = c{e,z2,-5/2,hl,^-i/2)ni-3,h2,e) 

+ ^ c{9, Z2, hi + 1/2, hl,^ - 1/2) {n{hi, h2, 9) - n{hi + 1, /i2, e)) , (6.10) 

fel . </ii<-4 

where, by definition, the function C{9, Z2, 171,171') is equal to 

^izi (h+ui N^)±iu!a {0-z)T/e O'l^l) 



C{9, Z2, 171, m') — e 



IZi 



)(f)i{9 - z) dzi. (6.11) 



Substituting (6.10) in (6.9) we obtain 
l^i^'^^l < Ci sup sup 



sup \Ci9,Z2,hi + l/2,hl,,-l/2)\ 

- <hi<-3 



x( ^ |7^(/il,-3,^) -7e(/ii + l,-3,^) 

+ \n{-3,h2,9)-n{-3,h2 + i,9) 



+ E E |7^(/ll,/l2,^)-7^(/ll + l,/l2,^) 

. <hi<-Ah? . <fe2<-4 

mm — ^ — mm — ^ — 

-n{hi,h2 + l,9) + n{hi + l,h2+l,9 
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By condition 14, all sums in (6.12) are bounded uniformly on 9 , since they do not exceed 
the variations of Ro(-, ^) on the set [ri — 4 — cr, ri + 4 + cr] x [r2 — 4 — cr, r2 + 4 + cr] . 
Hence, 

|/i'''^|<C2 sup sup sup |C(^, ;22,/ii + 1/2, /i^n- 1/2)1- (6.12) 

6le[-7r,7r]2 Z2e[-7r,7r] /i^.„<?ii<-3 

Lemma 6.2 ('see Lemma 3.7 from [3]) Let condition II (ii) hold. Then 

e^i'^-(^^'^^-'^)^/^C(e, Z2,m,m') 7r(signm - signm')0i(^i, ^2 - Z2), s ^ 0, (6.13) 

uniformly in 9 E [— 7r,7r]^ , Z2 G [— 7r,7r] and \m\, \m'\ > 2 . 

Since /ii + 1/2 < —2, ~ 1/2 ^ ~2 , the integral /J^'^-* vanishes as e — ^ by 

Lemma 6.2 and (6.12). Similarly, the remaining integrals in (6.5) with i = 1, 3 or j = 1, 3 , 
vanish as e — > , i.e., the series in I'^ over max^ \ hj\ > 3 vanish as £ — > . Proposition 6.1 
is proved. ■ 

Step (Hi) The next step in the proof is to prove the following asymptotics for the RHS of 
(6.4) as £ ^ : 



8 

[—n,7r]'' 



If = AZl^ j e-^^-('-^)+^'^<^'(^)^/^^(^)n,(^)Ro(rTVa;,(^)T,^)n,,(^) 



TT .-c /o _ .-t: /n - 



^Oi{z^)- (• • • ( / e-^('^+'^^^^)a(;.2) 

%Zd ^ IZ2 

— TT 

J5/2NeZi _ -i5/2NeZi . . . . 

^^z^{h+u,N,)±^u.^{e-z)r/e^(^^^^ _ _ _ _ )dZd)d9 + 0^(1) , (6.14) 

izi ) ) ) ) 

— TT 

where Ot-(1) — as e ^ for any r 7^ . Formula (6.14) was proved in [3, Lemma 3.8] 
for the case when d=\. This formula is based on the formula Yl,\mj-vj\<2F{/^3i^ei'^^ ~ 
^ivjN^zj (^^ib/2NeZj _ g-i5/2Are2j^ ^^^^ foUowiug inequality 

|Ro(£[r/£] - veN, - heN„ 9) - Ro(r ^ Vu;^{9)t, 9)\ < C{e + N,). 

This inequahty follows because \h\ < 2 and Ro(r, 6*) satisfies condition 14. The proof of [3, 
Lemma 3.8] admits extension to the case when d > 1 . 

Step (iv) Let us apply Lemma 6.2 to the inner integrals in the RHS of (6.14) over zi, . . . , 
and obtain 

If = i^!^ / e-'^-(^-^)e<'^-'('')±'^-(''))"/'^5(^)n^(^)Ro(r ^ Wuj^{9)t, 9)Ii^>{9) d9 + o(l)(6.15) 
8 J 

as £ ^ 0. Note that Ro^(r, •) e C(T'') , i,j = 0,1, by condition (2.15). Moreover, the 
identities ujcr{9) +cJo-'(^) = const+ or O0cr{9) — uo„i{9) = const_ with the const± 7^ are 
impossible by condition E5. Furthermore, the oscillatory integrals with uj^rif)) -^uo„i{9) ^ 
const± vanish as e — > by the Lebesgue-Riemann theorem. Hence, only the integrals with 
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uJaid) — uJcr'iO) = contribute to the integral (6.15), since Ua{9) +Ua'{9) = would imply 
LiJa{0) = = , which is impossible by E4. Thus, for cr 7^ cr' , = o{l) as £ — > . For 

a — a' , If — 0(1) and 

7- = j e-''^<^-p~^g{e)n^{e)tio{r + Vuj^{e)T, e)ii^,{e) de + o(i) as £ ^ o. 

This completes the proof of Lemma 3.6. ■ 



7 Appendix B: Proof of Lemma 3.7 

We first apply (3.18) and obtain 

S:,r/e = ^ E E Ql/eil + ^) E / s] - STuN,, y - x) sign([rjs] - y,)gi,,{p + yf . 

Let us write Kr,£,m — £[r/e] — smN^ . The Parseval equality yields 
E sign([ri/£] - yi)Ro(Kr,£,m, y - x)Qli^{p + yf 

= (27r)-'^ y F,_,, [sign([ri/£] - yi)Ro(«.,.,^, y - 2:)]f,_,,[^?^/^(p + y)^] dO' . (7.1) 



Note that 



1 \ 



i^.^.'[sign([n/5] - 2/0] = (27r)'^-^<5(^"0PV (^^^(^j ^'''^^^'^ 



where 9' — (^i, ^') , and PV stands for the Cauchy principal part. Hence, 

I'd \ 



SI.,. = -5(2.)--' E E g;,.(i+x)/ pv/ _,)^,) R.K......^-')& 



- —i^^r'-' EE/ e-^('+^>^^.^/.(^) 

me J xeImNs'Y2d 



J[ri/e] {d{-z)+ixiz 



I tgm - ^)/2) ^' ^"'^ '^^j de'dO. (7.2) 

We change variables: 9'i ^ (p = 9'i — z , and then denote 2; = 6*'^ . Therefore, 

S-r/s = EE/ e-^('-^--^')e-(^'-^)^,^/,(^)Ro(«,,., ^')^e(^') d^'d^, 

where 

h{e') := PV j Gr/M + ^, d'Tgio'i + ^, d^, e' = {e[, e'). 
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By formula (3.2), the matrix ^t(^) has the form 

gt(e) ^ cos Q{e)t I + sin n{e)tc {6), 

where / is the identity matrix and C{9) is introduced in (2.10). By Lemma 1.2 (iv), 



G:{0) = J2i<'OSu;,{0)t + smu;,{0)tC:{e))U,{0) 

a=l 

^ j2e^'^^('^-/'^-^^^^u,{e), eeT'xc,, (7.3) 



(T=l 



since cosoja{e)t = (e^'^"* + e-^-*)/2 and sina;^(^)i = {e^-^ - e-^-*)/{2i). Applying the 
partition of unity (3.5), (3.6), and formula (7.3) we rewrite leid) in the form 

4(^) = EEPV / g.{e,+v.ef- i^^^^^f±^n.(^,+^,^)d^. 



Lemma 7.1 (i) sup sup^^g l-^£(^)l < 

6»eT'',rieR 

(a) Let ViUJa{0) 7^ ±ri , for 6 e suppg^^ and for fixed ri e R. Then 

h{e) -2mJ2i2 9k{e)e^'^-^'^^/'f^,{e)^a{e) ^ O as e^+O, (7.4) 

k,± a=l 

for fixed r > and ri e R. ifere f^{0) is a matrix-valued function of the form 

/±(^) = sign(ri ± Viu;,(^)t)(/ ^ zC:(^))/2. 

Lemma 7.1 can be proved by using the technique of [8, Lemma 8.3] or of [1, Proposition 
A. 4 (i), (ii)]. The proof is based on the following well-known assertion 

lim I PV / '"'^^^'^^^"'^ dz - 7rze'^^(°)x(0) sign/'(0) | = 0, 

where x ^ C'S f & and /'(O) 7^ . Lemma 7.1 gives 

S-r/e - (2^)-''^^EEE E /^.(^)e-^('-^--^')e-(^'-^)^,^/,(^) 

a=l k,± me J xei"mJVt-rp2d 

xRo(Kr,,,m, e')e^'^^^''^^/' f^^{e')u,{e') de'de. (7.5) 



Comparing (3.21) and (7.5) we see that the problem of evaluating the limit value of (7.5) is 
solved by the similar way as in Lemma 3.6. ■ 
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8 Appendix C: Proof of Lemma 3.8 

By (3.20) we write 

Sl/e = E E E Sr/ei^ +x)R\Kr,e,m, [r/e]-X, [v j e\- y)Ql , + vf . 

where K,r,e,m = ^[^/^] ~ smN^ . Change variables y ^ z = y — x and denote the sum over 
m and x by ^eiz) , 

= E E G'r/.il+^)R'i^r,e,m,[r/e]-x,[r/e]-x-z)g'^/^{p + x + zf. (8.1) 
meJ xeimNs 

Therefore, 

Slr/e = E ^e{z). (8.2) 

The estimate (2.14) and definition (3.19) imply the same estimate for : 

\R%r,x,y)\<C{l + \x-y\r\ (8.3) 
Next, the Cauchy-Schwartz inequality yields 

E E Ksii+^)\Keip + ^ + ^f\ < Y.Ksii+^)\Ksip + ^ + ^f\ 

< ll^.%ll.'^<C(l + ||y-li.(T.)). (8.4) 
Hence, condition E6 and estimate (8.3) imply that |$£(-2)| < C(l + l-z])^'''. Since •j > d, 

y: m^)\<c<oc, (8.5) 

and the series in (8.2) converges uniformly in e (and also in r,r,l,p). Therefore, it suffices 
to prove that 

lim^J^) = for each z E Z'^. (8.6) 

Let us consider the series in (8.1). At first, note that by definitions 10 and (3.19), the 
function R^{r,x,y) depends on x — y , i.e., has the form R^{r,x,y) = R°(r, Xi, a; — y) , 
and R°(r, i/i, z) =0 for yi < . Hence, 

R°(r, [ri/e] — Xi, [ri/e] — Xi — Zi, z) = for xi > [ri/e] — Zi. 

Further, from condition (2.13) it follows that V5 > > such that for any yi > Ks 

|R'^(r, yi, yi — Zi,z)\ < S . Hence, V5 > 3Ms — max{Ks, Zi) > such that 

|E E Sr/ei^+^)^°i'^r,e,m,['ri/£]-Xl,[ri/e]-Xi-Zi,z)g^^^{p + X + z)^ 

me J x& I^N, 

xi < [ri/e] - Ms 

< ^E E \Q'r/eii+x)g'^/eiP + ^ + ^f\<CS, 

m&J X&ImNe 
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by estimate (8.4). It remains to prove that for fixed Ms > , 
E E G',/,{l+x)R''{Kr,e,m,[n/e]-xi,[n/e]-xi-zi,z)g'^/^{p + x + zf ^0, 8^0, (8.7) 

me J x&Am 

where Am = {x = {xi,x) : xi G {[ri/e] - Ms, [ri/e] - zi) fl ImiN,] x e IfaN,} ■ For enougfi 
small an e > 0, there is a G Z such that ([ri/e] — Ms, [ri/s] — Zi) C ImiN^ ■ Hence, 
|£:[ri/e] — em^^A^^I < C{eN^ + e) . Therefore, by condition 14, we have 

\R°{e[n/e] - em\N„ . . .) - i?°(0, . . .)| < Ce^"^ with some /3 G (0, 1). 

Hence, by the estimate (8.4), we can replace BP{e[ri/e\ — em\Ns, . . .) into BP{Q, . . .) in the 
series (8.7). 

Further, let us change Xi — >■ [ri/e] — Xi in (8.7). Therefore, Xi runs the finite number 
of points, Xi G {zi,Ms) . To derive (8.6) it suffices to prove that for every fixed z E Z"^ , 
xi eZ, r eW^ , and r 7^ , 

J2^\0,f-emN„Xi,Xi-Zi,z)p^ J2 Q^r/ei^ + Xr)apg',/AP + + z)^s ^ ^ (8-8) 

as £ — > , where a, P,j,S = 1, . . . ,d , Xy {\ti/£\ — Xi, x) . For simplicity of exposition, we 
omit indices i,j,k,l in (8.8) and assume that d — 1 . In this case, to prove (8.8) it suffices 
to show that 

g',/,{k + [ri/e] - x^)glj^{p + [n/e] - + z)^ ^ as £ ^ 0. 
Indeed, applying the Fourier transform, we have 

gi,^{h + [ri/£] - X,) = (27r)-^ j^^ e-'^'^+^^^'^^-^-^'^g,,,{e^)g{e^) dO,. (8.9) 

Using the decomposition (3.8) we rewrite (8.9) as 

e;,%(Zi + [ri/£]-xi) = (27r)-^ /^^.(^i)a±(^i)^(^i)e-^('i+['-^/^l-^^)^^±^'^'^(^^)^/^d^i. (8.10) 

The eigenvalues uJcf{0) and the matrices a^{6) are real-analytic functions inside the suppc/fc 
for every k . Moreover, conditions E4 and E6 imply that for fixed ri G and r 7^ , 
mcs {^^1 G : ViCJo-(^^i) = ±^i/t} = 0. Hence, the integrals in (8.10) vanish as e ^ by 
the Lebesgue-Riemann theorem. The proof of convergence (8.8) in the case when c/ > 1 is 
similar and based on condition 14. ■ 
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